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摘 要 : 为 分 析 研 究 生态 保护 与 修复 工程 实施 背景 下 甘肃 中 东部 植被 生长 季 NDVI 时 空 变化 及 其 对 
气候 因子 的 响应 ,采用 MODIS 数据 产品 MOD09Q1, 运 用 最 大 值 合成 法 、 趋 势 分 析 法 与 相关 分 析 法 ， 
对 该 区 域 2008—2016 年 NDVI 的 时 空 变化 特征 及 其 对 气候 因子 的 响应 进行 了 分 析 研 究 。 结 果 表 
明 :(1) 时 间 层 面 上 ,2008 一 2016 年 植被 生长 季 NDVI 呈现 出 增加 的 趋势 ,增加 速率 为 0.0014a， 


2012 年 增幅 最 大 ;空间 层面 上 ,2008 一 2016 年 研究 


区 内 的 低 山 丘陵 区 和 平原 区 的 植被 生长 季 NDVI 


增加 明显 ;部 分 中 起 伏 山地 植被 活动 减弱 。(2) 生长 季 NDVI 对 不 同 气 候 因 子 的 响应 程度 不 同 ,在 
P <0.05 的 显著 性 水 平 下 ,生长 季 NDVI 对 气温 的 变化 最 敏感 ,区 域内 大 部 分 面积 上 二 者 呈现 出 强 


相关 与 极 强 相关 ;降水 次 之 ,以 中 等 相关 与 强 相 关 为 主 ;太阳 辐射 最 小 ,以 弱 相 关 
研究 结果 可 为 针对 区 域 差 异性 并 有 所 侧重 的 生态 工程 的 实施 提供 相应 的 参考 与 依据 ,从 而 加 强 区 


其 至 不 相关 为 主 。 


域 生 态 安全 ,进而 促进 国家 生态 安全 屏障 的 加 快 建设 。 


关键 ia: 


近年 来 ,改善 生态 环境 以 缓解 现 阶段 社会 主要 
矛盾 备 受 关注 '""。 作 为 陆地 生态 系统 敏感 成 员 的 
植被 能 够 通过 其 最 直观 的 表现 较 好 的 表征 区 域 生态 
环境 状态 2 ,其 覆盖 变化 能 够 反映 一 定 区 域内 自 
然 、 人 工 和 经 济 生态 系统 中 各 因素 对 植被 的 影响 程 
度 -5 。 采 用 人 工 治理 修复 与 自然 恢复 相 结合 生态 
保护 与 修复 工程 耦合 了 自然 和 社会 经 济 生态 系 统 ， 
使 得 区 域 的 整体 生态 功能 逐步 恢复 和 提升 。 目 前 关 
于 生态 工程 方面 的 研究 主要 包括 生态 工程 治理 成 
效 “-" 及 其 效益 评价 * -中 方面 和 工程 实施 现状 及 
问题 分 析 "等 方面 。 气 候 是 决定 地 球 上 植被 覆盖 
状况 的 最 主要 因素 ,作为 自然 因素 的 气候 因子 与 植 
被 间 相互 关系 的 研究 已 成 为 全 球 变化 研究 的 重要 内 
容 呈 ,而 定量 指标 归 一 化 植被 指数 NDVI 与 植被 覆 
盖 之 间 呈 正 相关 关系 握 - 59 ,日 被 广泛 用 于 植被 变化 
和 植被 一 气候 关系 的 研究 。 因 此 ,可 运用 植被 生长 
Æ NDVI 来 探析 生态 工程 实施 背景 下 区 域 植被 覆盖 
变化 状况 及 其 对 气候 因子 的 响应 情况 ,从 而 实现 自 


Q@” 收 稿 日 期 : 2019-04-05; 修订 日 期 : 2019 -07 -10 


NDVI; 趋势 分 析 法 ; 气候 因子 ; 甘肃 中 东部 


然 与 人 为 双重 因素 对 生态 环境 影响 的 初步 探究 。 目 
前 关于 植被 生长 季 NDVI 应 用 方面 的 研究 主要 集中 
在 区 域 和 项 目 区 上 ,具体 体现 为 区 域 尺度 上 植被 覆 
盖 时 空 变化 及 驱动 因素 分 析 方 面 “” "| 和 项 目 尺度 
上 单一 研究 对 象 的 变化 方面 “" 。 基 于 生态 工程 
背景 对 植被 生长 季 NDVI 时 空 变化 及 其 对 气候 因子 
响应 方面 的 研究 鲜 有 涉及 。 

考虑 到 黄土 高 原 地 区 相对 恶劣 的 自然 条 件 和 脆 
弱 的 生态 系统 ,区 域内 严重 的 水 土 流 失 及 其 治理 成 
BE SATE” ,植被 状况 作为 生态 治理 成 效 最 为 直 
观 的 体现 亦 给 予 了 高 度 关注 。 甘 肃 中 东部 地 区 脆弱 
的 生态 环境 一 直 以 来 都 是 生态 工程 治理 的 观点 所 
在 ,近年 来 对 于 该 区 域 植被 覆盖 变化 方面 的 研究 鲜 
有 涉及 。 基 于 此 ,本 研究 以 2008 一 2016 年 为 研究 时 
段 ,采用 MOD09Q1 数据 ,基于 ENVI Classic 5.3 + 
IDL 8.5 与 ArcGIS 10.2 -F AA MRT 工具 对 研究 时 
段 内 的 WD 了 打数 据 进 行 相应 的 处 理 , 来 探究 多 种 生 
态 工程 人 工 治理 修复 与 自然 修复 对 区 域 植被 覆盖 状 
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况 的 影响 ,并 进一步 探究 其 对 气候 因子 的 响应 ,从 而 
分 析 科 学 评价 研究 区 生态 环境 质量 提升 程度 ,实现 
多 年 来 黄土 高 原 地 区 生态 工程 治理 成 效 的 定性 分 
析 , 为 后 期 区 域 生态 保护 与 修复 政策 的 制定 提供 相 
应 的 参考 ,实现 对 区 域 生 态 环 境 保 护 与 修复 的 进 一 
步 指导 。 


1 数据 来 源 与 研究 方法 


1.1 研究 区 概况 

本 研究 中 的 甘肃 中 东部 地 区 为 典型 的 黄土 丘陵 
ARK ,具体 的 地 理 位 置 如 图 1 所 示 ,涉及 甘肃 省 陇 
中 、 东 地 区 和 陇南 山地 ,共计 4 个 市 区 。 其 中 , 庆 阳 
市 与 平凉 市 属于 生态 环境 极为 脆弱 . 沟 坝 地 和 坡 旱 
地 广泛 分 布 的 陇 东 地 区 ;中 部 偏 南 的 定西 市 属于 黄 
土 高 原 西 部 边缘 地 带 和 西 秦 岭 末 端的 陇 中 地 区 ;天 
水 市 则 由 于 北 秦 岭 与 祁连山 的 共同 作用 ,属于 区 域 
内 多 分 布 向 斜 构 造 和 断 陷 谷 地 的 陇南 山地 。 无 论 是 
植被 还 是 降水 分 布 ,整个 研究 区 均 因 黄土 高 原 独 特 
的 气候 而 呈现 自 东 南 至 西北 地 带 性 分 异 的 规律 。 自 
1999 年 退耕 还 林 还 草 生态 工程 的 实施 ,区 域内 的 植 
被 状况 有 所 好 转 , 至 2012 年 东部 百 万 亩 土地 整治 重 
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大 工程 项 目 以 及 其 后 实施 的 山水 林 田 湖 草 生态 保护 
与 修复 工程 , 均 使 得 区 域内 生态 环境 状况 明显 改善 ， 
人 民生 活 质量 再 创新 高 , 旨 在 加 强 区 域 生 态 环境 改 
善 .保护 与 管 护 ,以 及 更 进一步 提高 人 民生 活水 平 的 
重大 工程 相继 开展 。 主 观 上 生态 环境 得 到 极 大 改 
善 ,但 从 植被 角度 来 客观 分 析 生 态 工 程 治 理 成 效 方 
面 的 研究 相对 缺乏 。 
1.2 数据 来 源 与 预 处 理 

本 研究 中 2008 一 2016 年 的 MODIS ( Moderate- 
Resolution Imaging Spectroradiometer ) 数据 产品 
MOD09Q1 来 源 于 美国 LP DAAC (Land Process Dis- 
tributed Active Archive Center , USA) ,主要 采用 USGS 
EROS 数据 中 心 所 开发 的 MRT (MODIS Reprojection 
Tool) 工 具 来 完成 对 影像 子 波段 (第 1、2 波段 ) 的 格 
式 与 投影 转换 以 及 重 采 样 等 处 理 。DEM 数据 来 源 
于 美国 地 质 勘 探 局 ,是 ASTER GDEM V2 30 m 分 辨 
率 数 字 高 程 。 降 水 .气温 与 太阳 辐射 数据 来 源 于 寒 
区 旱 区 科学 数据 中 心 ,SWAT 模型 中 国 大 气 同化 驱 
动 数 据 集 (CMADS V1.0), MODO9Q1 数据 产品 空 
间 分 辩 率 为 250 m, 时间 分 辨 率 为 8 d。 本 研究 主 
要 基于 ENVIClassic5.3 +IDL8.5 及 其 插件 “小 能 
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图 1 地 理 位 置 区 位 图 


Fig.1 Location map 
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TAA V3.0”、ArcGIS 10.2 与 MRT 工具 来 对 研究 
区 的 原始 影像 进行 处 理 , 其 中 采用 最 大 值 合成 法 获 
取 了 研究 区 植被 生长 季 的 月 度 NDVI 数据 。 由 于 影 
像 拍摄 过 程 中 会 受到 云 太阳 高 度 角 等 的 影响 ,影像 
数据 会 产生 异常 值 , 需 要 对 其 进行 滤波 处 理 , 本 研究 
中 采用 的 是 S-G( Savtzky-Golay ) 滤波 ?| 法。 

1.3 研究 方法 

1.3.1 最 大 值 合成 法 (MVC) 本 研究 中 的 最 大 值 
合成 法 MVC(Maximum Value Composites ) (22) 主要 用 
于 植被 生长 季 (4 ~ 10 月 ) NDVI 月 度数 据 的 获取 。 
该 方法 能 够 在 一 定 程度 上 减少 原始 影像 拍摄 过 程 中 
云 、 太 阳 高 度 角 等 的 影响 ,其 计算 公式 如 下 所 示 : 
MaxNDVI, = Max ( NDVI, ,NDVI,,NDVI,,NDVI,) (1) 


式 中 :MaxNDVI; 表示 月 度 NDVI 最 大 值 ;i 为 月 份 ， 
取 值 为 4~10;NDVI ~ NDVI, 分 别 表示 相应 月 份 内 
原始 影像 预 处 理 后 的 NDVI 值 。 

1.3.2 Savtzky-Golay 滤波 法 ”Savizky-Golay 滤波 
法 即 Savizky-Golay 滤波 拟 合 法 , 它 通 过 确定 合理 的 
滤波 参数 ,对 研究 区 植被 生长 季 NDVI 数据 作 加 权 
多 项 式 拟 合 , 以 求 得 最 小 均 方 根 误差 ,这 样 远离 大 多 
数 点 的 一 些 边缘 点 便 不 会 参与 拟 合 一 。 故 拟 合 时 
过 于 偏离 正常 生长 趋势 线 的 噪声 部 分 会 被 丢弃 。 其 
基本 公式 如 下 所 示 : 


m 
$, C x Ypa 
_ i=-m 


y* 
N 


(2) 
SUP Y, Y 分 别 表示 原 始 NDVI 数值 和 平滑 后 的 
NDVI 数值 ;C; 为 多 项 式 拟 合 系数 ,表示 从 滤波 器 首 
部 开始 的 第 i 个 NDVI (ERILE; N 为 滤波 器 的 长 
度 ,其 大 小 等 于 2 m+1。 

1.3.3 趋势 分 析 法 ”趋势 分 析 法 ”的 原理 在 于 运 
用 一 元 线性 回归 分 析 的 方法 对 区 域内 时 序 植被 生长 
Æ NDVI 的 变化 趋势 和 其 强度 进行 分 析 。 其 中 , 回 
归 方 程 的 斜率 slope 代表 研究 时 段 内 每 个 栅 格 点 植 
被 指数 的 变化 趋势 。 若 slope > 0 ,代表 该 栅 格 植被 
状况 有 所 改善 , 且 slope 值 越 大 改善 效果 越 明 显 ; 反 
之 亦 然 ; 若 slope =0, 则 表示 植被 状况 无 变化 。 


nx 3 (i x Myon, ) - D iÈ Moon, 
` a. (3) 


nx Dë- (Di 
式 中 ;slope 是 趋势 线 的 斜率 ;n 是 时 间 序 列 的 长 度 ， 


slope = 


本 研究 中 为 9 a(2008 一 2016 年 ) ;i 是 第 i 年 ;NDVI 
是 第 i 年 生长 季 NDVI 值 (4 ~10 月 )。 


2 结果 与 分 析 

2.1 甘肃 中 东部 2008 一 2016 年 植被 生长 季 NDVI 
时 空 变 化 分 析 

2.1.1 甘肃 中 东部 2008 一 2016 年 植被 生长 季 ND- 


WI 年 际 变 化 分 析 ”本 研究 首先 基于 8 d 合成 的 
MOD09Q1 数据 的 第 1、2 波段 计算 得 到 了 相应 时 长 
研究 区 的 NDV 数据 ,然后 运用 最 大 值 合成 法 计算 
得 到 研究 区 植被 生长 季 (4 ~ 10 A) NDVI 月 度数 据 ， 
最 后 得 到 研究 区 2008 一 2016 年 年 际 NDVI 均值 及 
其 变化 情况 如 图 2 所 示 。 

由 图 2 可 知 ,总 体 上 植被 生长 季 NDVI 年 际 变 
化 呈现 出 增加 趋 热 , 且 以 0.001 4. a ! 的 速率 增加 ， 
2012 年 是 一 个 明显 的 突变 点 ,增幅 最 大 。2008 一 
2011 年 ,生长 季 年 际 NDVI 均值 变化 相对 稳定 ,保持 
在 0.191 9 ~0.199 4 之 间 。 这 主要 是 由 于 人 研究 区 所 
属 区 域 为 典型 的 黄土 压 陵 沟 蜜 区 ,生态 环境 相对 脆 
弱 ,多 年 来 由 于 退耕 还 林 还 草 等 各 种 生态 工程 的 洽 
理 , 生 态 环境 状 况 有 所 改善 ,人 民生 活水 平 有 所 提 
升 ,其 对 于 生态 环境 的 认 知 程度 亦 得 到 有 效 改 观 ,最 
为 显著 的 是 2008 一 2011 年 甘肃 省 广泛 开展 的 退耕 
还 林 还 草 工程 成 效 的 保存 率 和 管 护 率 均 达 到 99% 
以 上 ,种 植 结构 中 饲 草 类 作物 的 种 植 比例 提高 了 
10% ,区 域内 生态 治理 成 效 显著 。 

2012 年 作为 一 个 明显 的 突变 点 , 可 以 看 出 
2011 一 2012 年 ,生长 季 年 际 NDVI 变化 最 大 ,达到 了 
0.030 6, 2011 年 是 研究 时 段 内 实施 退耕 还 林 还 草 
面积 最 大 的 一 年 ,2008 年 开展 的 大 规模 土地 整治 特 
别 注重 研究 区 内 时 坡地 的 整理 ,于 滞后 期 2 a 后 土 
地 生产 力 显著 提高 ,导致 2012 年 NDVI 变化 最 大 。 
2012 年 后 ,生长 季 NDVI 保持 基本 稳定 。 这 主要 是 
由 于 2013 年 研究 区 内 开展 了 甘肃 东部 百 万 亩 土地 
整治 重大 项 目 工 程 ,主要 针对 研究 区 内 的 黄土 梁 闭 
与 黄土 卉 等 水 土 流 失 较 为 严重 的 地 区 进行 了 相应 的 
治理 ,于 2014 和 2015 年 再 次 开展 的 新 一 轮 退耕 还 
林 还 草 , 尤 其 注重 特色 林果业 的 发 展 ,注重 生态 效益 
与 经 济 效益 兼 优 ,并 优先 安排 生态 环境 脆弱 、 生 态 地 
位 重要 的 区 域 ,同时 加 强 退 耕地 保护 与 退耕 经 营 权 
适度 流转 , 均 使 得 区 域内 的 土地 生产 能 力 不 仅 得 到 
提高 ,农民 的 收入 也 显著 增加 。2015 年 后 各 地 充分 
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2 甘肃 中 东部 2008 一 2016 年 植被 生长 季 NDVI 年 际 变化 


Fig.2 Variation of annual NDVI of vegetation growing season in mid and eastern Gansu Province from 2008 to 2016 


认识 到 山水 林 田 湖 是 一 个 生命 共同 体 ,对 山水 林 田 
湖 草 进行 统一 保护 与 修复 必须 遵循 自然 规律 ,使 得 
区 域内 NDVI 变化 相对 稳定 并 有 所 提高 。 
2.1.2 甘肃 中 东部 2008 一 2016 年 植被 生长 季 ND- 
Viz team 本 研究 基于 ENVI5.3 平台 ,采用 
时 序数 据 趋势 分 析 模 型 ,对 研究 区 9 a 的 NDVI 空间 
变化 情况 进行 了 相应 分 析 ,具体 结果 如 图 3 所 示 。 
图 3a 为 2008 一 2016 年 NDVI 空间 变化 ,变化 量 
介 于 -0.290 8 ~0. 2699 之 间 。 陇 东 地 区 的 NDVI 
的 增加 状况 明显 较 好 , 主要 以 平凉 与 庆 阳 两 市 的 交 
界 处 最 为 明显 ,而 平凉 东部 以 及 庆 阳 南部 NDVI 的 
变化 相对 较 小 ;隶属 于 陇南 山地 的 天 水 市 的 NDVI 
变化 呈现 出 自 东 南 向 西北 递减 的 趋势 ,以 东部 NDVI 
增加 最 为 明显 ; 陇 中 地 区 的 NDVI 变化 则 较 小 ,尤其 
是 定西 中 部 地 区 呈现 出 NDVI 负 向 增加 情况 ,但 总 
体 上 NDVI 仍 是 增加 的 。 以 上 变化 均 说 明 研 究 区 内 
所 实施 的 各 种 生态 保护 与 修复 工程 过 多 注重 各 自治 
理工 程 的 成 效 , 忽 视 了 各 工程 相互 之 间 的 联系 ,导致 
区 域 整体 上 生态 治理 成 效 欠 佳 。 近 年 来 以 “山水 林 
田 湖 是 一 个 生命 共同 体 ” 的 倡议 为 指引 ,区 域内 开 


台 注重 生态 系统 内 部 之 间 的 有 机 联系 ,生态 环境 的 
全 面 改 善 与 修复 有 待 进一步 发 展 。 
图 3b 为 2008—2016 年 植被 指数 变化 趋势 的 笠 
Z (slope) , 介 于 -0.042 3 ~0.038 4 之 间 , 说 明 植 被 
变化 情况 相对 较 大 。 依 据 罗 敏 等 所 的 划分 依据 ,说 
明 区 域内 植被 变化 相对 稳定 ,基本 无 退化 现象 。 
slope 值 比较 高 的 区 域 主要 分 布 在 以 平凉 与 庆 阳 两 
市 交界 处 为 代表 的 低 山 丘陵 区 和 平原 地 区 ,这 些 地 
区 植被 生长 季 NDVI 明显 上 升 ,植被 覆盖 程度 增加 ， 
植被 活动 明显 增强 。slope 值 较 低 的 区 域 主要 分 布 
在 以 平凉 东部 .定西 与 天 水 两 市 交界 处 .定西 中 西部 
地 区 以 及 庆 阳 南部 的 黄土 坊 和 中 起 伏 山地 ,这 些 地 
区 植被 生长 季 NDVI 有 所 下 降 , 植 被 覆盖 程度 下 降 ， 
植被 活动 相对 减弱 。 植 被 活动 减弱 可 主要 归 因 于 区 
域内 相对 较 差 的 自然 条 件 , 年 降水 量 明 显 低 于 蒸发 
量 , 加 之 生态 工程 的 实施 破坏 了 原 有 的 平衡 系统 ,区 
域 植被 耗 水 与 蒸 散发 加 剧 ,导致 土壤 出 现 干 层 , 出 现 
部 分 植被 死亡 ,导致 植被 NDVI 指数 相应 有 所 下 降 ， 
植被 活动 减弱 。 人 研究 区 内 并 无 slope 为 0 的 区 域 。 
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图 3 甘肃 中 东部 2008 一 2016 年 NDVI 空间 变化 及 变化 趋势 图 


Spatial and trend variation of NDVI in mid and eastern Gansu Province from 2008 to 2016 


Fig. 3 


赵 鸿雁 等 :甘肃 中 东部 植被 生长 季 NDVI 时 空 变化 及 其 对 气候 因子 的 响应 


2.2 甘肃 中 东部 2008—2016 年 植被 生长 季 NDVI 
对 气候 因子 的 响应 分 析 

由 于 植被 覆盖 变化 对 气候 变化 的 响应 在 不 同 区 
域 上 差异 较 大 ,研究 区 属于 干旱 半 干 旱地 区 , 故 气 
温 \ 降 水 与 太阳 辐射 会 直接 或 间接 地 对 区 域 植被 覆 
盖 状 况 造 成 影响 。 气 温 上 升 增强 了 植物 内 部 酶 的 活 
性 ,加 速 了 化 学 反应 速率 ,使 得 植被 春秋 季 物 候 期 
有 所 提前 或 推迟 ,延长 了 植被 生长 季 ,植被 覆盖 状况 
有 所 改善 ;气温 过 高 会 使 植物 的 蒸腾 作用 加 快 , 导 致 
植被 耗 水 与 蒸 散发 量 均 增 加 ,土壤 水 分 加 快 流失 , 植 
被 活动 因 水 分 限制 而 相应 减弱 。 而 降雨 量 会 直接 影 
响 土壤 含水 量 ,同样 会 加 速 植物 内 部 化 学 反应 速 
度 ' 站 。 太 阳 辐 射 较 好 的 阳 坡 由 于 降雨 量 较 少 会 容 
易 导 致 更 多 的 薰 散 发 ,使 植被 覆盖 降低 , 阴 坡 则 由 于 
光照 和 太阳 辐射 较 弱 而 容易 形成 干 冷 的 天 气 ,植被 
活动 明显 较 弱 ,其 至 不 利于 植被 的 生长 '*]。 

因此 ,为 探究 生态 保护 与 修复 工程 实施 情况 下 
气候 因子 (月 累积 降水 /mm .月 平均 气温 /YX 与 月 平 
均 太阳 辐射 /MJ © m”) 对 植被 生长 季 NDVI 的 影 
响 ,本 研究 先 从 时 间 层 面 上 分 析 了 年 均 NDVI 与 气 
候 因子 的 变化 情况 如 图 4 所 示 。 

由 图 4 可知, 生长 季 NDVI 与 降水 的 变化 情况 
除 2010 年 .2015 年 与 2016 年 外 ,其 余年 份 变化 趋 
势 基本 一 致 ,表现 出 一 定 程 度 上 的 同步 。 生 长 季 
NDVI 与 气温 变化 趋势 除 2012 年 .2016 年 外 ,其 余 
年 份 变化 趋势 基本 一 致 。 生 长 季 NDVI 与 太阳 辐射 


0.23 p 600 
0.22 + 500 
z 0.21 + 400 和 
人 020 | 300 5 
0.19 | 200 se 
0.18 f 一 一 NDJT ”一 口 一 降水 100 


0.17 0 


Co (=A Le = N oD = ite) Ke] 
© O 一 一 一 一 一 一 一 
=] So (= © © (=) (=) So © 
N N N N N N N N N 
年 份 
0.23 上 760 
0.22 上 740 E 
z 0.21 上 720 3 
z 0.20 | = 
= 0.19 | og = 
0.18 - —o— NDVI 一 上 太阳 辐射 680 & 
K 
0.17 L L L 1 上 1 L 660 
oo © SS oe N oD < 10 Ko] 
© oS -| m= Do | Bown Dem | re Do | 
© © © So © © © © © 
N N N N N N N N N 


变化 除 2011 年 .2014 年 与 2015 年 外 ,其 余年 份 变 
化 趋势 基本 吻合 。 总 体 来 说 , 较 降 水 与 太阳 辐射 , 气 
温 与 生长 季 NDVI 的 变化 趋势 吻合 程度 更 高 。 

鉴于 年 际 气 候 因子 与 生长 季 NDVI 呈现 出 基本 
一 致 的 变化 趋势 ,为 进一步 探讨 二 者 之 间 相 关 关 系 
的 密切 程度 ,本 研究 基于 像 元 尺度 ,运用 小 能 工具 箱 
v3.0 中 的 时 序数 据 相 关 分 析 模 型 来 探究 2008 一 
2016 年 空间 上 气候 因子 与 植被 生长 季 NDVI 间 相 关 
关系 的 密切 程度 。 时 序数 据 相 关 性 分 析 模 型 是 用 来 
分 析 两 组 多 波段 数据 之 间 的 相关 性 ,输出 两 组 多 波 
段 数据 的 相关 系数 波段 。 本 研究 以 2008 一 2016 年 
生长 季 NDVI 月 度数 据 为 输入 数据 ,输出 2008 一 
2016 年 生长 季 NDVI 与 各 气候 因子 之 间 的 相关 系 
数 。 根 据 相 关系 数 的 等 级 划分 ,将 其 分 为 基本 不 相 
关 (1rl 和 0.2)、 弱 相关 (0.2 < 1r1 和 0.4) .中 等 相关 
(0.4 < Irl <0. 6) . 强 相关 (0.6 < Irl <0. 8) 、 极 强 相 
关 (0.8 < Irl <1)5 级 ,并 统计 得 到 2008 一 2016 年 
不 同 相 关 ( 忆 <0. 05) 等 级 下 生长 季 NDVI 与 降水 、 气 
温和 太阳 辐射 之 间 的 占 比 情况 如 图 5 所 示 。 

由 图 5 可 知 ,不 同 相关 等 级 下 ,降水 与 生长 季 
NDVI 之 间 在 区 域内 大 部 分 面积 上 显示 为 中 等 程度 
相关 与 强 相 关 , 弱 相关 与 基本 不 相关 次 之 , 极 强 相 关 
Ae, A SAK NDVI 之 间 在 区 域内 大 部 
分 面积 上 显示 为 强 相 关 与 极 强 相关 ,中 等 程度 相关 
面积 量 占 比 次 之 ,基本 不 相关 与 弱 相 关 面 积 占 比较 
> ,太阳 辐射 与 生长 季 NDVI 之 间 在 区 域内 大 部 分 
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图 4 甘肃 中 东部 2008 一 2016 年 生长 季 NDVT 与 
气候 因子 年 际 变化 
Fig.4 Slope map and annual variation of different 
slopes’ NDVI in mid and eastern Gansu Province 
from 2008 to 2017 
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图 5 甘肃 中 东部 2008—2016 年 生长 季 NDVI 与 气候 因子 间 的 相关 等 级 占 比 


Fig.5 Proportion of related grades between NDVI and climatic factors in mid and eastern Gansu Province from 2008 to 2016 


面积 上 显示 不 相关 或 者 弱 相 关 , 少 部 分 区 域 显示 中 
等 程度 相关 , 极 少 部 分 区 域 显示 出 强 相 关 与 极 强 相 
关 。 因 此 ,总 体 来 说 ,气温 对 植被 的 生长 状况 最 敏 
感 ,降水 次 之 ,太阳 辐射 最 小 。 

综 上 ,2010 年 .2011 年 .2012 年 .2014 年 .2015 
年 与 2016 年 生长 季 NDVI 变化 与 气候 因子 变化 不 
一 致 。2010 年 区 域内 生长 季 NDVI 与 降水 之 间 的 相 
关系 数 以 强 相关 与 极 强 相关 为 主 ,但 该 年 份 生长 季 
NDVI 较 2009 年 有 所 下 降 , 降 水 是 上 升 的 ,而 气温 是 
显著 下 降 的 ,太阳 辐射 是 增加 的 ,这 进一步 说 明生 长 
Æ NDVI 对 气温 的 敏感 性 ,气温 的 下 降 导致 区 域内 
以 农作物 为 主 的 部 分 植被 长 势 欠 佳 ,从 而 导致 生长 
Æ NDVI 下 降 。2011 年 气温 与 降水 均 呈 现下 降 趋 
势 , 尽 管 太 阳 辐 射 上 升 不 少 ,但 生长 季 NDVI 还 是 呈 
现 出 下 降 的 趋势 ,这 也 说 明 降 水 与 气温 对 生长 季 
NDVI 影响 较 大 ,太阳 辐射 影响 较 小 。2012 年 降水 
和 和 气温 与 生长 季 NDVI 的 相关 性 较 强 ,在 各 气候 因 
子 均 有 所 增加 的 情况 下 ,区 域内 生长 季 NDVI 明显 
增加 。2014 年 与 2015 年 各 气候 因子 均 呈 直线 下 降 
的 趋势 ,区 域内 大 部 分 面积 的 生长 季 NDVI 与 气温 
和 降水 呈现 出 强 相 关 与 极 强 相关 ,可 见 该 阶段 NDVI 
下 降 与 气温 和 降水 关系 密切 。2016 年 降水 有 所 回 
升 ,但 是 年 均 气温 升 高 了 1 % ,气温 与 生长 季 NDVI 
之 间 的 关系 以 极 强 相关 为 主 ,气温 升 高 导致 区 域内 
部 植被 欠 住 ,导致 区 域内 生长 季 NDVI 下降。 总体 


来 说 ,尽管 研究 区 属于 干旱 半 干 旱 区 ,降水 并 不 是 影 
响 区 域内 植被 生长 状况 的 最 主要 因素 ,这 与 区 域 地 
处 黄土 丘陵 沟 密 区 有 极 大 的 相关 性 ,一 系列 生态 工 
程 对 水 土 流失 等 的 治理 使 得 降水 与 植被 生长 状况 不 
再 呈现 出 负 相关 关系 ,并 逐渐 转向 正 相 关 。 植 被 生 
长 对 气温 的 变化 仍 是 最 敏感 的 。 


3 讨论 


赵 安 周 . 石 玉 琼 和 秦 国 玲 等 研究 得 到 大 型 生态 
工程 对 NDVI 均 会 产生 正面 影响 "|。 甘 肃 中 东 
部 地 区 土 层 厚 .地 IIR ARIE ICE RACH , 生 
态 环境 相对 脆弱 。1999 年 以 来 生态 工程 通过 人 工 
种 植 ,组 建新 的 植被 系统 ,打破 了 原 有 的 平衡 系统 ， 
植被 耗 水 大 于 降水 补给 ,强化 了 土壤 干燥 化 程度 。 
针对 这 一 现象 ,从 土壤 水 分 的 角度 来 说 ,在 生态 工程 
实施 初期 ,土壤 条 件 较 好 ,植被 树木 较 多 , 枯 枝 落叶 
和 腐殖质 增多 ,土壤 表面 性 质 得 到 改善 ,土壤 蕾 水 量 
增强 ,表层 土壤 水 分 有 所 增加 2 。 生 态 工程 实施 
6 ~12 a 间 存 在 大 幅 波 动 ,12 a 之 后 灌木 就 会 成 为 工 
程 区 的 优势 种 ,群落 结构 趋 于 稳定 , 动 植物 残 体 逐渐 
增多 ,改善 了 土壤 环境 ,使 得 土壤 水 分 逐渐 增多 , 植 
被 覆盖 相应 增加 '*”1。12 a 之 后 随 着 土壤 深度 的 增 
加 ,土壤 含水 量 减 少 的 越 多 ,导致 土壤 干 层 出 现 。 故 
对 于 植被 覆盖 度 合 理 阔 值 的 确定 以 及 下 垫 面 特征 如 
土壤 水 分 对 植被 覆盖 的 影响 有 待 进一步 研究 。 


赵 鸿 雁 等 :甘肃 中 东部 植被 生长 季 NDVI 时空 


本 研究 得 到 区 域内 2008 一 2016 年 植被 生长 季 
NDVI 呈现 出 明显 增加 的 趋势 ,这 与 誉 丹 丹 等 ”| 、 曹 
博 等 ”的 研究 结果 一 致 , 且 ND 增加 区 域 集中 在 
陇 中 和 陇 东 地 区 。 区 域内 植被 生长 季 NDVI 对 气候 
因子 中 的 气温 最 敏感 ,二 者 之 间 的 相关 性 在 95% 的 
著 性 水 平 下 大 部 分 以 强 相关 和 极 强 相 关 为 主 ; 降 
水 次 之 ,以 中 等 相关 与 强 相关 为 主 ;太阳 辐射 的 影响 
则 为 最 小 ,以 弱 相 关 甚 至 不 相关 为 主 。 这 与 曹 永 强 
等 "1 曹 博 等 .的 研究 结果 相似 ,鉴于 区 域 间 的 差 

异性 ,相关 性 程度 也 有 所 不 同 。 故 气候 因子 对 生长 
Æ NDVI 影响 程度 的 定量 研究 有 待 进一步 探讨 。 


Si 
i: 


4 结论 


(1) 甘肃 中 东部 地 区 2008 一 2016 年 年 际 NDVI 
在 时 间 层 面 上 呈现 为 增加 的 趋势 ,增加 速率 为 
0.001 4 . a ,2012 年 增幅 最 大 ;在 空间 层面 上 ,其 


空间 变化 图 与 slope 图 均 显 示人 研究 区 内 的 低 山 丘陵 
区 和 平原 区 NDVI 增加 最 为 明显 ,植被 活动 明显 增 
强 ;而 在 研究 区 内 的 部 分 中 起 伏 山地 ,植被 活动 出 现 
减弱 的 现象 ,但 后 期 正确 认识 到 了 多 种 生态 工程 之 
间 的 有 机 联系 ,生态 治理 后 的 管 护 力度 有 所 增加 , 植 
被 活动 相应 增强 ,生态 环境 质量 明显 提升 。 

(2) 生 长 季 NDVI 对 不 同 气候 因子 的 响应 程度 
不 同 , 其 中 气温 变化 与 生长 季 NDVI 变化 关系 最 密 
切 ,在 区 域内 大 部 分 面积 上 呈现 出 强 相关 与 极 强 相 
关 , 降 水 次 之 ,以 中 等 相关 与 强 相 关 为 主 ,太阳 辐射 

的 影响 则 为 最 小 ,以 弱 相 关 其 至 不 相关 为 主 。 本 研 
究 对 研究 区 植被 生长 季 NDVI 的 时 空 变化 及 其 与 气 
候 因 子 变化 之 间 的 相关 性 所 进行 的 研究 可 为 针对 区 
域 差 异性 并 有 所 侧重 的 生态 工程 的 实施 提供 相应 的 
参考 与 依据 ,从 而 加 强 区 域 生态 安全 ,进而 促进 国家 
生态 安全 屏障 的 加 快 建设 。 
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Spatiotemporal variation of NDVI in vegetation growing season 
and its responses to climatic factors in mid and eastern Gansu 
Province from 2008 to 2016 


ZHAO Hong-yan', CHEN Ying”, ZHOU Yi’, PEI Ting-ting’ , 
XIE Bao-peng , WANG Xiao-qian' 
(1 College of Resources and Environment , Gansu Agricultural University , Lanzhou 730070 ,Gansu, China; 
2 College of Management , Gansu Agricultural University , Lanzhou 730070 , Gansu, China) 


Abstract: As a sensitive member of terrestrial ecosystem , vegetation can characterize the influence extent of natu- 
ral and human factors on the regional ecological environment. In view of the scarcity of good ecological environment 
in the world, ecological engineering was used as a means to protect and restore the ecological environment. The 
quantitative index normalized difference vegetation index (NDVI) was widely used to study vegetation change and 
the relationship between vegetation and climate. Therefore ,in order to analyze the NDVI variation in the vegetation 
growing season in mid and eastern Gansu Province , northwest China under the background of implementing the Eco- 
logical Engineering for Ecological Protection and Restoration, MODO09Q1 data was used to analyze the spatiotemporal 
variation of NDVI in the vegetation growing season and its responses to climatic factors by applying maximum value 
composite (MVC) ,trend and correlation analysis approach from 2008 to 2016. The results showed as follows: (1) 
On the time level ,the NDVI of the vegetation growing season showed an increasing trend from 2008 to 2016, with an 


! and the largest increase occurred in 2012. On the spatial level ,the NDVI of the vege- 


increasing rate of 0.001 4 a` 
tation growing season in the study area, the low hilly areas and plains, increased significantly from 2008 to 2016. 
Furthermore , the vegetation activity in the rolling hills in some areas is weakened. (2) The response degree of ND- 
VI to different climatic factors is different in vegetation growing season. At the significant level of P <0. 05, NDVI 
variation in the growing season is the most sensitive to the change of temperature. Most of the areas in the region 
show strong and extremely strong correlation,and the precipitation is the second sensitive factor with moderate and 
strong correlation,and the solar radiation is the least, with weakly correlation or even without correlation. The re- 
search can provide the corresponding reference and basis for the implementation of the ecological engineering based 
on regional differences and the particular purposes ,so as to strengthen the regional ecological security and accelerate 
the construction of the national ecological security barrier. 


Key words: NDVI; trend analysis approach; climatic factors ;mid and eastern Gansu Province 


